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on  a  few  populations  early  and  possibly  make  the  final  decision  before  the  censoring 
time.  An  example  is  provided  in  the  final  part  to  illustrate  the  use  of  the  early  selection 
rule  for  Weibull  populations. 
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§  1  Introduction 


Consider  designing  and  analyzing  an  experiment  for  comparing  k  populations  tti, 
7r2,  •  •  • ,  TTfc-  Suppose  that  m  items  are  taken  from  each  population  and  observations  can 
be  obtained  from  those  items  in  time  order,  as  for  example,  in  a  life-testing  experi¬ 
ment.  It  is  often  desirable  to  terminate  the  test  from  a  population  as  soon  as  there  is 
enough  statistical  evidence  that  it  is  not  the  best  population,  and  then  this  population 
is  eliminated  from  further  consideration. 

Assume  that  the  random  observations  from  population  ttj  have  a  density  function 
f{x\9i)  of  the  form 

f{x\9)  =  c{0)exp{P{6)Q{x)}h{x),  xeX,  (1.1) 

where  X  is  the  support  of  f{x\9)  and  h{x)  >  0  for  x  £  X.  Let  ft  be  the  parameter  space 
for  each  9i. 

Let  0[i]  <  9[2]  <  •  •  •  <  %]  denote  the  ordered  values  of  the  parameters  9i,92,-  "  i^k- 
It  is  assumed  that  the  exact  pairing  between  the  ordered  and  the  unordered  parameters 
is  unknown.  The  population  associated  with  the  largest  value  9[k]  is  considered  as  the 
best  population.  The  readers  are  referred  to  Gupta  and  Panchapakesan  (1979)  for  a 
comprehensive  understanding  of  selection  and  ranking  procedures. 

The  function  f{x\9)  is  assumed  to  have  (nondecreasing)  monotone  likelihood  ratio 
with  respect  to  x.  This  assumption  is  equivalent  to 

im)  -  m)){Qi^2)  -  Q{xi))  >  0  (1.2) 

for  any  02  >  9i,  0i,02  E  X2  >  xi,  Xi,X2  E  X.  Without  loss  of  generality,  we  assume 
that  Q{x)  is  a  nondecreasing  function  of  x  and  /3(0)  is  a  nondecreasing  function  of  9. 

Many  exponential  family  distributions,  such  as  Chi-square,  Exponential,  Gamma  (o:, 
P)  with  one  of  the  two  parameters  known.  Log-normal  {p,  a^)  with  a  known,  Weihullij, 
P)  with  one  of  the  two  parameters  known,  have  property  (1.2).  So  our  results  here  can 
be  applied  to  them. 

In  an  application  situation  of  industrial  life-testing  experiment,  m  items  from  each 
of  the  k  population  tti,  •  •  •,7rfc  are  independently  put  on  test  at  the  outset  and  are  not 
replaced  on  failure.  Due  to  the  time  restriction,  the  experiment  terminates  at  a  pre¬ 
specified  time  T.  The  failure  time  of  an  item  is  observable  only  if  it  fails  before  time 
T.  Otherwise  the  item  is  said  to  be  censored  at  time  T.  This  type  of  time  censoring  is 
known  as  the  type-I  censoring.  The  type-I  censoring  scheme  has  received  much  attention 
in  the  statistical  literature,  see  Spurrier  and  Wei  (1980),  and  others.  The  ranking  and 
selection  procedures  based  on  censored  data  for  the  exponential  distribution  have  been 
considered,  for  example,  in  Berger  and  Kim  (1985),  Gupta  and  Liang  (1993),  and  Huang 
and  Huang  (1980). 

In  this  paper,  we  derive  a  Bayes  selection  rule  for  exponential  family  distributions 
based  on  type-I  censored  data.  A  monotone  property  of  this  rule  is  discussed  and  an 
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early  selection  rule  is  proposed.  Through  this  early  selection  rule,  one  can  terminate  the 
experiment  on  a  few  populations  early  and  possibly  make  the  final  decision  before  the 
given  time  T.  The  approach  used  here  is  similar  to  that  of  Gupta  and  Liang(1993). 


§2  A  Bayes  Selection  Rule 


Let  T  be  the  censoring  time.  Let  Xij,  1  <  j  <  m  be  the  type-I  censored  data  of  the 
m  items  taken  from  population  tt*.  We  only  observe  •min{Xij,T).  Let  Ni  =  I[Xij<T] 
be  the  number  of  uncensored  data  up  to  time  T. 

Let  Yii  <  Yi2  <  ■  •  ■  <  Yii^^  be  the  ordered  values  of  the  Ni  observed  data  given  Ni. 
Then  (V^i,  1^2)  •  •  * ,  Yini^,Ni)  have  a  joint  probability  density 


fiVunilOi) 


m\ 


_  n  hiVij)  (2.1) 

{m  ni)i 


ml 


[m  —  rii)] 


where  yi  =  {yn,  yi2,  •  •  • ,  yin),  0  <  n  <  m,  <  2/^2  <  •  •  •  <  ym  <  T  and 

QiVij)  +  {m-  ni)Q{T).  (2.2) 

j=i 

_  Let  Af  be  the  sample  space  generated  by  A'’  =  (ui,n2)  •  •  •  ,«jfc)  and  conditioned  on 

N  =  n  =  (ni,  n2,  •  •  • ,  njt),  let  be  the  sample  space  generated  by  T  =  (yi,  y2,  •  •  • ,  y^)- 
Let  0  =  (0i,02,  -  ■  ■ , 6k)  and  Q  =  {d\6i  e  fi,  1  <  i  <  A:}  be  the  parameter  space.  Let  A 
be  the  action  space.  Action  i  corresponds  to  the  selection  of  population  TTj  as  the  best 
population.  For  a  given  6  E  fl  and  an  action  i,  the  associated  loss  function  is  defined  by 


L*{e,i)^L{e^k]-6i)  (2.3) 

where  L{x)  is  a  nonnegative  and  nondecreasing  function  of  x,  a:  >  0,  such  that  L(0)  =  0. 

Let  g{6)  =  n)=i  gj{6j)  be  the  prior  density  over  the  parameter  space  Q.  It  is  assumed 
that  L(6ik])g(6)d6j<  00.  . 

A  selection  rule  S  =  (61,62,'  ‘ '  ,^k)  is  defined  to  be  a  measurable  mapping  from  the 
sample  space  (Af,  (yn)ne//)  [^>  ^  ^  ^j=i  ^)  =  1  for 

all  y  E  yn,n  e  Af.  The  value  of  6i{y,  n)  is  the  probability  of  selecting  population  Wj  as 
the  best  p(y)ulation  based  on  the  observation  (y,n). 

Let  R(6,g)  denote  the  Bayes  risk  associated  with  the  selection  rule  6.  Then  by 
Fubini’s  Theorem  we  have 

k 

R{^^g)=Y^  [  [  L{6[k]- 6i)f{y,n\9)g(d)dMy  (2.4) 
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where  f{y,n\$)  =  Ili^i /(j/i,  Now  let 

MVi:  rii)  =  f{yi,  ni\9i)gi{ei)d6i,  f{y,  h)  =  fj  fi{yi,  Ui), 

and  g{d\y,n)  =  l[gi{ei\yi,ni). 

Then  (2.4)  becomes 

k 

9)  =  ^  f  f~  “  ^i)90\y,  n)def{y,  h)dy. 

For  each  {y,n),  define 

Ai{y,  ^)  =  J~  ^(%]  “  ^i)9{^y,  n)de,  i  =  1, 2,  •  •  • ,  fc, 


and  let 


^(y,«)  =  {*|Ai(y,n)  =  min  Aj{y,n)}. 

1<3<K 

Then  a  uniformly  randomized  Bayes  rule  is  So  =  (^gii  •  •  •  >  ^Gk),  where 

f  \My^^)\~^  if  ieA{y,n) 

I  0  otherwise. 


(2.5) 


(2.6) 


(2.7) 


§3  A  Monotonicity  Property  of  6c 


For  each  fixed  {yuUi),  gi{Bi\yuni)  =  0  if  and  only  if  ^i(^i)  =  0.  Then  gi{di\y*,ni)  and 
gi^ilyuni)  have  the  common  support.  Let  Di  be  their  common  support.  Consider  the 
likelihood  ratio  defined  on  Di,  by 


riiOi\y*,n*,yi,ni) 


9i{0i\yt,n*i) 

9i{0i\yi,ni) 


(3.1) 


A  simple  calculation  shows  that  for  some  nonnegative  function  W 


ri{0i\y*,n*,yi,ni) 

roo 

=  exp{m)lQ(x)  -  Q{T)]}k{x)dzY'‘'-"’\ 
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from  which  we  get  the  following  lemma. 

Lemma  3.1  Let  ri{6i\y*,n*,yi,ni)  be  defined  by  (3.1).  Then 

(a)  forn*i  =  n^y*  >  yi,  ri{9i\yf,n*i,yi,ni)  is  a  nondecreasing  function  of  Oi  in  Di  and 

(b)  foryl  =  >  ni,  ri{9i\y*,n\,yi,ni)  is  a  nonincreasing  function  of  9i  in  Di. 

The  following  lemma  is  used  in  the  proof  of  lemma  3.3. 

Lemma  3.2  If  g{9)  and  h{9)  are  probability  density  functions  such  that  g{9)/h{9)  is 
nondecreasing  function  of  9  in  f2,  then  for  any  nonincreasing  function  f{9)  of  9  in  D, 

f  f{9)h{9)d9  >  f  f{9)g{9)d9. 


Lemma  3.3  Let  Ai{y,h)  be  defined  in  (2.5).  For  each  i(l  <  *  <  k),  Ai{y,n)  is 
nonincreasing  in  yi  and  also  in  nj,j  ^  i  when  all  the  other  variables  are  kept  fixed,  and 
nondecreasing  in  ni  and  also  in  yj,j  ^  i,  when  all  the  other  variables  are  kept  fixed. 


Proof.  We  only  prove  that  Afiy,  h)  is  nonincreasing  in  yi  when  all  the  other  variables 
are  kept  fixed.  The  other  parts  can  be  proved  in  a  similar  way. 

Define 

9^  =  (^1 ,  •  ■  ■  1  )  ^i+1 )  ■  ■  ■  ) 

D*  =  {F  :9j  eO., j  =  1,2,- ■•,k,ji^i} 

y  = 

y  ~  iy\ )  ■  ’ ' )  2/i— 1 )  yi  >  yi+x >  ■  ■  * )  yfc)  • 


Then 


Since  for  each  fixed  0*  and  n,  L{9[k]  -  9i)  is  nonincreasing  in  9i  and  by  Lemma  3.1, 
ri{9i\yi,n*,yi,ni)  is  a  nondecreasing  function  of  9i  for  y*  >  yi.  So  Lemma  3.2  implies 
that 

J  9i)gi(9i\yi,ni)d9i  ^  j^L{6^k^  —  9i)gfi9i\yi,ni)d9i 

and  hence  Ai{y,n)  >  Ai{y*,h). 


Now,  from  Lemma  3.3,  we  obtain  a  monotone  property  for  Sq  in  the  following  theo¬ 
rem. 


Theorem  3.4  For  each  i  —  1,2,  Saijj.h)  is  nondecreasing  in  yi  and  nonin¬ 

creasing  in  ni ,  when  all  the  other  variables  are  kept  fixed. 
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Proof.  We  only  prove  that  8Gi{y,n)  is_nondecreasing  in  yi  when  all  other  variables 
are  kept  fixed.  The  monotone  property  of  5g  in  rij  can  be  proved  in  a  similar  way. 

Use  the  notation  in  Lemma  3.3.  Assume  y*  >  yi. 

If  i  ^  A{y,h),  then  5Gi{y,n)  —  0.  Since  SGi{y*,n)  is  nonnegative,  SGi{y*,n)  > 
8Gi{y,  n). 

If  i  €  A(y,  n),  then  Ai{y,n)  <  minj^^j  Aj(y,  n).  Using  Lemma  3.3, 

Ai(y*,n)  <  Ai{y,n)  <  min  Aj(y,n)  <  minAj(y*,n). 

And  hence  i  G  A{y*,n). 

To  get  5Gi{y*,n)  >  5Gi{y,n),  we  still  need  to  show 

A{y*,n)  C  A{y,n).  (3.2) 

For  each  h  G  A{y*,n),  Ah{y*,n)  —  Ai{y*,n).  Using  Lemma  3.3, 

^h{y,n)  <  Ah{y*,n)  =  Ai(y*,n)  <  Ai{y,n)  =  min  Aj{y,n) 

and  hence  h  G  A{y,n).  So  (3.2)  is  proved  and  SGi{y*,n)  >  SGi{y,fi). 


§4  An  Early  Selection  Rule 

In  this  section,  we  consider  the  following  linear  loss  function:  —  $i)  =  6^^]  —  9i, 

the  difference  between  the  parameters  of  the  best  and  the  selected  populations.  Thus 
the  set  A{y,  n)  given  by  (2.6)  turns  out  to  be: 

=  {^1  /  di9i{Gi\yuni)d9i  =  max^.y^  ^i9j{^j\yj',nj)d9j).  (4.1) 

Similar  to  the  proof  of  Lemma  3.3,  we  can  prove  the  following  result. 

Lemma  4.1  For  each  fixed  i,  E[9i\yi,ni]  is  increasing  in  yi  and  decreasing  in  Ui. 

Now,  we  will  use  Lemma  4.1  to  derive  an  early  selection  rule. 

At  time  t,  0  <  t  <  T,  let  Nfit)  denote  the  number  of  uncensored  data  from  population 
TTj  upto  t.  That  is,  Ni{t)  =  :  1  <  i  <  m.Xij  <  t}.  Also,  let  <  yi2  <  •  •  •  < 

^iNiit)  denote  observed  uncensored  data  given  Nfit).  At  time  t,  we  can  make  early 
decision  as  follows: 

Declare  population  tt;  as  a  non-best  population  and  exclude  it  from  further  experiment 
if  there  exists  some  population  tt/,  such  that 

Nh{t)  <  m  and  E{9h\yh{t),m\  >  E[9i\yi{t,T),Ni{t)\  (4.2a) 

or 

Nh{t)  =  m  and  E[9h\yh{t),m]  >  E[9i\yi{t,T),  Nfit)]  (4.26) 
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where 


Nkit) 

yh{t)  =  Y.  QiVhj)  +  {'m-  Nh{t))Q{t)  (4.3a) 

i=i 

and 

Ni{t) 

y,(t,T)  =  Y,  Q(v.i)  +  (™  -  rnWin  (4.36) 

Let  S{t)  denote  the  indices  of  the  contending  populations  for  the  best  at  time  t.  That 
is, 

S{t)  =  :  Nh{t)  <  (=)m  and  El0ilyi(t,T),  Ni(t)]  >  (>)El&klyh(t),m],  h  ^  i}.  (4.4) 

The  following  lemma  shows  that  for  any  t,  Q  <  t  <T,  S{t)  is  not  empty. 

Lemma  4.2  For  any  Q  <  t  <T,  the  set  S{t)  defined  by  (4-4)  is  not  empty. 

Proof.  Let 


S'{t)  =  {i :  E[ei\yi{t),  Ni{t)]  =  m&x  E[eh\yh{t),  Nh{t)]}. 

l<n<K 

Then  S'{t)  is  not  empty.  We  prove  that  S'{t)  is  a  subset  of  S{t).  For  i  E  S'{t)  and  any 

h^i, 

if  Nh{t)  <  m,  then 


E[eMt,T),Ni{t)]  >  E[9i\yi{t),Ni{t)]  >  E[eh\yn{t),  N^it)]  >  E[en\yh{t),m]- 
if  Nh{t)  —  m,  then 

E[ei\yi{t,T),Ni{t)]  >  E[ei\yi{t),Ni{t)]  >  E[e^\yn{t),  N^it)]  >  E[9h\yh{t),m]. 

In  either  situation,  we  see  that  i  G  S{t).  Hence  S'{t)  C  S{t). 

Now,  the  experiment  terminates  as  soon  as  there  is  a  time  t,  0  <  t  <  T,  such  that 
\S{t)\  —  1  and  in  this  case,  we  select  the  population  with  its  index  in  S{t)  as  the  best 
population.  Otherwise,  the  experiment  goes  on  until  time  T.  Let 

S{T)  =  {i  ;  E[9i\yi,  Ni]  =  max  E[9j\yj,  Nj]},  (4.5) 

where  S{T~),  which  is  not  empty  by  Lemma  4.2,  denotes  the  set  of  the  indices  of  those 
populations  having  not  been  eliminated  before  time  T.  Then,  a  uniformly  randomized 
selection  is  made  from  S(T). 

^From  the  above  description,  we  see  that  the  early  selection  rule  can  possibly  make 
a  final  selection  earlier  than  the  termination  time  T.  Denote  this  early  selection  rule  by 
—  (^Gii  ■  ■  ■  ’  ^hk)-  Then,  we  have  the  following  theorem. 
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Theorem  4.3  Under  the  loss  function  L{9),  5*q^  =  6Gi{y,n)  for  all  1  <  i  <  k,  y  E  yn 
and  h  G  Af,  where  dciiy,  h)  is  defined  by  (4-1)  and  (2.7). 

Let  ti  =  inf{t :  |5(i)l  =  1, 0  <  t  <  T}  A  T,  where  aAb  —  min{a,  b).  Then  Theorem 
4.3  is  equivalent  to  the  following  theorem. 

Theorem  4.4  S{ti)  =  A{y,n}  for  all  {y,h). 

Proof.  Case  1.  If  <i  <  T,  then  \S{ti)\  =  1.  Without  loss  of  generality,  we  let  TTjt  be 
the  population  with  index  in  the  set  S{ti).  Since  A{y,h)  contains  at  least  one  element,  it 
suffices  to  show  that  i  ^  A{y,h)  for  all  i  k.  Since  i  ^  S{ti),  it  means  that  population 

TTj  is  eliminated  at  some  prior  time,  say  to,  That  is,  at  time  to,  for  some  tt/,,  either 

Nh{to)<m  and  ElOhlyhito),^]  >  E[9i\yi{to,T),  Ni{to)]  (4.6a) 

or 

Nh{to)='m  and  E[9h\yhito),m]  >  E[9i\yi{to,T),  Ni{to)].  (4.66) 

Now,  note  that  Ni{t)  is  a  nondecreasing  function  of  t  G  (0,  T]  and  Ni(t)  <  m.  Also,  by 
(4.3a)  and  (4.3b),  yh{t)  is  nondecreasing  in  t  and  yi{t,T)  is  nonincreasing  in  t.  Especially, 
we  have 


Nh  =  Hh{T)  <  m,  Ni{t)  <  Ni{T)  =  W,  yi{to,T)  >  yt 


and 


yh 


>  yh{to)  if  Nk{to)  <  m, 
=  yh{to)  i{Nh{to)=m. 


Thus,  when  Nh{to)  =  m,  then  Nh  =  m.  Then  by  Lemma  4.1  and  (4.6b), 


E[9h\yh, 


=  E[9h\yh{to),m] 

>  E[9i\yi{to,T),  Ni{to)] 

>  E[9i\yi,Ni]. 


When  Nh{to)  <  rn,  then  yh  >  yhik)  and  Nh  <  m,  Therefore,  by  Lemma  4.1  and 
(4.6a), 


E[9h\yh,Nh]  =  E[9h\yh{to),'m] 

>  E[9i\yi{to,T),Ni{to)] 

>  E[9i\yi,Ni]. 

In  either  situation,  we  see  that  i  ^  A{y,n). 

Case  2.  If  =  T,  we  need  to  prove  that 

(a)  i  ^  S{T)  implies  i  ^  A{y,n),  and  (b)  i  €  S{T)  implies  i  e  A{y,n). 
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We  prove  (a)  first.  Suppose  i  ^  5(T).  Then,  TTj  is  eliminated  at  a  time  <T  hy 
some  other  tt/j. 

If  to  <  T,  this  reduces  to  the  situation  discussed  in  Case  1. 

If  to  =  T,  then  by  (4.5),  E[9h\yh,Nh]  >  E[9i\yi,Ni].  Therefore,  by  the  definition  of 
i  ^  A{y,n). 

For  (b),  we  have  firstly  A{y,n)  C  S{T)  C  S{T~)  by  (a)  and  definition  of  S{T)  and 
S{T-).  If  ieS{T-), 


E[9i\yi,Ni 


=  .^^^_^E[9j\yj,Nj] 

>  m^  E[9j\yj,Nj] 

jeA{y,n) 


This  means  i  £  A{y,n).  The  proof  now  is  completed. 


§  5  An  Example 

We  use  the  simulated  data  to  illustrate  how  the  early  selection  rule  works.  Suppose 
that  we  have  five  populations  TTj,  t  =  1,  2,  3,  4,  5.  The  lifetime  of  the  population  tti 
follows  a  Weibull  distribution  with  density 

^  >  0- 

The  unknown  parameters  9i,  9^  are  simulated  independently  from  17(0, 1).  That  is, 

9i,  •••,  9^  are  independent  and  identically  distributed  with  17(0, 1).  Ten  observations 
are  simulated  independently  from  each  population.  The  data  are  listed  in  the  following 
table. 


TTi  7r2  TTs  7r4  IT5 


0.10 

0.33 

0.50 

0.17 

0.50 

0.13 

0.13 

0.28 

0.36 

0.10 


0.06 

0.09 

0.05 

0.04 

0.02 

0.02 

0.07 

0.10 

0.07 

0.07 


0.57 

0.56 

0.61 

0.75 

0.68 

0.22 

0.47 

0.45 

0.69 

1.01 


0.67 

0.90 

1.02 

0.64 

0.74 

0.90 

0.94 

0.86 

0.59 

0.90 


0.80 

1.19 

0.70 

0.36 

0.88 

0.52 

0.71 

0.77 

0.52 

0.24 
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We  want  to  select  a  population  with  the  largest  mean  lifetime.  Since  the  mean  lifetime 
of  the  population  ttj  is  proportional  to  6i,  what  we  need  to  do  is  to  find  a  population 
with  the  parameter  9\k]-  Suppose  that  the  type-I  censoring  scheme  is  planned  before  the 
life-testing  experiment  and  the  censoring  time  is  set  to  be  T  =  1.  Therefore  we  obtain 
the  following  table. 


TTi  7^2  TTs 

E[ei\yi,ni]  0.14  0.03  0.31  0.34  0.83 


According  to  the  selection  rule  Sg,  at  the  end  of  the  experiment,  we  select  tts  as  the 
best  population. 

However,  if  the  early  selection  rule  Sq  is  applied,  we  can  make  the  selection  before 
T  =  1  and  end  the  experiment  earlier.  According  to  the  selection  rule  Sq,  at  time  t, 
0  <  t  <  T  =  1,  exclude  the  population  tt^  as  a  non-best  population  and  remove  it  from 
further  experiment  if  there  exists  some  population  tt/,  such  that 

Nh{t)  <  m  and  E[Gh\yh{t)-,m]  >  E[9i\yi{t,T),Ni{t)] 


or 

Nh{t)  =  m  and  E[9h\yh{t),'nA  >  E[9i\yi{t,T),Ni{t)]. 

According  to  this  rule,  at  ti  =  0.88,  all  the  populations  tti,  7r2,  tts  and  7r4  are  removed 
from  the  experiment  and  the  population  tts  is  selected  as  the  best.  So  the  experiment 
can  be  ended  at  t\  =  0.88  and  the  time  saved  is  0.12  or  12%. 
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comments  and  suggestions. 
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